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Abstract. In order to understand a crucial role of orbital degree of freedom in the 
magnetic structure of recently synthesized neptunium compounds NpTGa5 (T=Fe, 
Co, and Ni), we propose to discuss the magnetic phase of an effective two-orbital 
model, which has been constructed based on a j-j coupling scheme to explain the 
magnetic structure of uranium compounds UTGas. By analyzing the model with 
the use of numerical technique such as exact diagonalization, we obtain the phase 
diagram including several kinds of magnetic states. An orbital-based scenario is 
discussed to understand the change in the magnetic structure among C-, A-, and G- 
type antiferromagnetic phases, experimentally observed in NpFeGas, NpCoGas, and 
NpNiGa 5 . 
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1. Introduction 

Recently it has been widely recognized that orbital degree of freedom plays an essential 
role for the understanding of novel magnetism in transition metal oxides pQ. A 
typical material is manganese oxide, exhibiting remarkable colossal magneto-resistance 
phenomena |2j- Due to competition and interplay among spin, charge, and orbital 
degrees of freedom, rich phase diagram has been revealed [I], but a recent trend is 
to unveil further new phases both from experimental and theoretical investigations. In 
fact, even in undoped RMn03 (R=rare earth lanthanide ions), a novel antiferromagnetic 
(AF) phase called the "E-type" spin structure has been reported as the ground state for 
R=Ho [313]. Here we follow the definitions of spin structure in Ref. [5]. The origin of 
the E-AF phase has been clarified theoretically jE] based on a band-insulator scenario 
in the e g -orbital systems El E] ■ It should be noted that the ground state of undoped 
manganites was just considered to be well understood, since for most R-ions, the A-type 
AF insulating phase appears with the C-type ordering of the (3x 2 — r 2 )- and (3y 2 —r 2 )- 
orbitals [TUl [TH H21 CS1 [13 [lEl HZl CHI CHI 1201 12H 123 Moreover, also for 
the half-doped manganite Lao.sCao.sMnOs, a charge-ordered ferromagnetic (FM) phase 
has been found in experiments EH] , as predicted theoretically [23 ESI ■ These facts 
clearly indicate the importance of both experimental and theoretical efforts to unveil 
new phases in manganites in addition to the explanation of the complex phases already 
observed. Such efforts have also been made to find new phases in other transition metal 
oxides, for instance, ruthenates [2HI ED] and nickelates [HI] • 

A trend to seek for new magnetic as well as superconducting phases has been also 
found in the /-electron system, which is another type of spin-charge-orbital complex. 
Among many kinds of /-electron materials, in recent years, /-electron compounds with 
HoCoGa 5 -type tetragonal crystal structure [see Fig. 1(a)], frequently referred to as 
"115", have been intensively investigated both in experimental and theoretical research 
fields of condensed matter physics. Such vigorous activities are certainly motivated 
by high superconducting transition temperature T c observed in some 115 compounds. 
Especially, amazingly high value of T C =18.5K has been reported in PuCoGa 5 [321EII1E3] 
and the mechanism has been discussed theoretically [S3 EE]- Among 115 compounds, 
interesting magnetic properties have been reported for UTGas, where T is a transition 
metal ion 03 [M EOl EII E2 ESI Ell E3 EHI EZl EE] In particular, neutron scattering 
experiments have revealed that UNiGas exhibits the G-type AF phase, while UPdGas 
and UPtGa 5 have the A-type AF state [HJ EH]- Note that G-type indicates a three- 
dimensional Neel state, while A-type denotes a layered AF structure in which spins align 
ferromagnetically in the ah plane and AF along the c axis [Sj. It is quite interesting 
that the magnetic structure is different for U-115 compounds which differ only by the 
substitution of transition metal ions. 

Quite recently, 115 compounds including neptunium have been skillfully synthesized 
and several kinds of physical quantities have been successfully measured [IHl E01 EH E2 
EH1 EU E3 EH] ■ Especially, the de Haas- van Alphen (dHvA) effect has been observed in 
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(a) (b) NpFeGa 5 (c) NpCoGa 5 (d) NpNiGa 5 
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Figure 1. (a) Crystal structure of NpTGas. Schematic views of magnetic structures 
composed of magnetic moments of Np ions for (b) NpFeGas, (c) NpCoGas, and (d) 
NpNiGas. For NpFeGas, magnetic moments at Fe sites are also depicted. 



NpNiGas [SO] , which is the first observation of dHvA signal in transuranium compounds. 
For NpCoGas, the dHvA oscillations have been also detected and a couple of cylindrical 
Fermi surfaces are found [52] • For NpFeGas, the magnetic moment at Fe site has been 
suggested in neutron scattering experiments faSj and it has been also detected by 57 Fe 
Mossbauer spectroscopy [S3]. Interestingly enough, the magnetic structure of Np-115 
compounds also depends sensitively on transition metal ion [SHIES]: C-AF for NpFeGas, 
A-AF for NpCoGa 5 , and G-AF for NpNiGa 5 , as shown in Figs. l(b)-(d). Note that for 
NpNiGas, the G-AF structure is composed of canted Np moments and the peak in the 
neutron scattering intensity grows after the FM transition occurs [SHI EE] ■ In any case, 
it is characteristic of U-115 and Np-115 compounds that the magnetic properties are 
sensitive to the choice of transition metal ions. 

The appearance of several kinds of AF states reminds us of the magnetic phase 
diagram of manganites and thus, we envisage a scenario to understand the complex 
magnetic structure of actinide compounds based on an orbital degenerate model similar 
to that of manganites. However, one must pay close attention to the meanings of "spin" 
and "orbital" in /-electron systems. Since they are tightly coupled with each other 
through a strong spin-orbit interaction, distinguishing them is not straightforward in 
comparison with ci-electron systems. This point can create serious problems when we 
attempt to understand microscopic aspects of magnetism and superconductivity in f- 
electron compounds. Thus, it is necessary to carefully define the terms "orbital" and 
"spin" for / electrons in a microscopic discussion of magnetism and superconductivity 
in actinide compounds. 

In order to overcome such problems, we have proposed to employ a j-j coupling 
scheme to discuss /-electron systems [HS]- Here we stress the advantages of the j-j 
coupling scheme. First, it is quite convenient for the inclusion of many-body effects 
using standard quantum-field theoretical techniques, since individual /-electron states 
are clearly defined. In contrast, in the LS coupling scheme we cannot use such standard 
techniques, since Wick's theorem does not hold. Second we can, in principle, include 
the effects of valence fluctuations. In some uranium compounds, the valence of the 
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uranium ion is neither definitely U 3+ nor U 4+ , indicating that the /-electron number 
takes a value between 2 and 3. In the j-j coupling scheme this is simply regarded as 
the average number of / electron per uranium ion. 

In this paper, we attempt to explain the complex magnetic structure of Np-115 
based on the effective two-orbital model, which has been constructed from the j-j 
coupling scheme in order to understand the magnetic properties of U-115 compounds 
|57j . We depict the Fermi surfaces of the kinetic term of 5/ electron based on the j-j 
coupling scheme in comparison with those obtained from the dHvA experiments. The 
agreement seems to be fairly well, but the present itinerant picture should be critically 
discussed. Then, we apply the exact diagonalization technique to the model in a small 
2x2x2 cube to obtain a clue to understand the complex magnetic structure of NpTGas. 
The phase diagrams are found to include G-, A-, and C-type AF states, consistent with 
experimental observations in Np-115. 

The organization of this paper is as follows. In Sec. 2, we will introduce an effective 
Hamiltonian, called the Tg model, for actinide 115 systems. In Sec. 3, we show the 
calculated results of the model Hamiltonian. The Fermi-surface structure is discussed 
in comparison with dHvA experimental results. We show numerical results for the 
magnetic structure obtained by using exact diagonalization technique. Finally, in Sec. 4, 
future developments are discussed and the paper is summarized. Throughout the paper, 
we use units such that h—k-Q—1. 

2. Model Hamiltonian 

In this section, we show an effective model based on the j-j coupling scheme for Np-115 
compounds, which is the same as the model for U-115 jHZj- We emphasize that the 
model with orbital degree of freedom is applicable to actinide 115 materials in common. 

2.1. Local f -electron state 

In order to construct a microscopic Hamiltonian for /-electron systems, it is necessary 
to include simultaneously the itinerant nature of / electrons as well as strong electron 
correlation and the effect of crystalline electric field (CEF). For the purpose, we have 
proposed to use the so-called j-j coupling scheme As shown in Fig. 2, we include 
the spin-orbit coupling so as to define the state labeled by the total angular momentum. 
For / orbitals with £=3, we immediately obtain an octet with j=7/2(=3+l/2) and a 
sextet with j=5/2(=3— 1/2), which are well separated by the spin-orbit interaction. In 
general, the spin-orbital coupling is as large as leV in actinide elements and thus, in 
the j-j coupling scheme, we consider only the j—5/2 sextet. 

In actual compounds, due to the electrostatic potential from ligand ions surrounding 
the actinide ion, the six-fold degeneracy of j—5/2 is lifted, as shown in Fig. 2. Note 
that in the j-j coupling scheme, it is enough to define the one-electron potential, which 
is deduced from the CEF level scheme of corresponding /^electron system. First 
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Figure 2. Schematic view for local /-electron state in the j-j coupling scheme. The 
j=5/2 sextet is split into several Kramers doublets due to cubic or tetragonal crystalline 
electric field effect. 

we consider the one /-electron state in the AuCu 3 -type cubic crystal structure, since 
this is the mother structure of HoCoGas-type tetragonal compound. The effects of 
tetragonality will be discussed later. For the case of cubic symmetry, we identify two 
eigen energies as — 240-8° for the T 7 doublet and I2OS4 for the T 8 quartet, where B® is 
a cubic CEF parameter. For the following discussion, it is quite convenient to introduce 
"pseudospin" to distinguish the degenerate states of the Kramers doublet and "orbital" 
to label the different Kramers doublets in r 8 quartet. In Fig. 3(a), we show the shape 
of the wavefuction of two r 8 orbitals, r 8 and r 8 . 

In order to proceed with the discussion, we need to know which is lower, or r 8 , 
in the one /-electron picture. For some crystal structures, it is possible to determine the 
level scheme from intuitive discussions on /-electron wavefunctions and the positions 
of ligand ions. However, this is not the case for the AuCu3-type crystal structure. For 
this case we invoke certain experimental results on Celn 3 , a typical AuCu3-type Ce- 
based compound, where T 7 and T 8 have been reported as ground and excited states, 
respectively, with an energy difference of 12meV [58- Thus, we take T 7 to be lower for 
the present considerations, as shown in Fig. 3(b). 

In the j-j coupling scheme for UGa3 and NpGa3, it is necessary to accommodate 
three or four electrons in the one-electron energy states T 7 and r 8 . We immediately 
notice that there are two possibilities, i.e., low- and high-spin states, depending on the 
Hund's rule interaction and the splitting between the T 7 and T 8 levels. As discussed in 
Ref. pSj, the effective Hund's rule interaction becomes small in the j-j coupling scheme 
and thus, the low-spin state should be realized, as shown in Figs. 3(c) and (d). We 
emphasize that the low-spin state is consistent with the LS coupling scheme. In the 
electron configuration shown in Figs. 3(c) and (d), the T 7 level is fully occupied to form 




Figure 3. (a) Two orbitals, a and b, in the r§ model. Level scheme for (b) Celn3, 
(c) UGa3, and (d) NpGa3 based on the j-j coupling scheme. Here we assume trivalent 
actinide ions as 

tj3+ ( 5 j3) and Np 3+ (5/4) It should be noted that up and down 
arrows denote pseudospins to distinguish the states in the Kramers doublet. Note also 
that for NpGa3, a couple of electrons in Tg orbitals form a local triplet, leading to IV 



a singlet. If this level is located well below the r 8 , the occupying electrons will not 
contribute to the magnetic properties. Thus, we ignore the electrons in the following 
discussion. As for details to validate the suppression of T 7 level, readers should refer 
Ref. We also mention another theoretical effort based on the j-j coupling scheme, 
in which the dual nature, itinerant and localized, of 5/ electron has been emphasized 
[531 EH EH- 

As shown in Fig. 1(a), since 115 compounds have the tetragonal structure, we need 
to include the effect of tetragonality. In the tetragonal structure, the quartet r 8 is 
further split into two doublets, Tq and I\, as shown in Fig. 2. Namely, there appear 
three Kramers doublets with one r 6 and two V-?. Note that two states are obtained 
from the mixture of and Tg in the cubic symmetry. However, for simplicity, we ignore 
the change of wavefunctions from cubic to tetragonal systems, since we believe that the 
existence of orbital degree of freedom is a key issue to understand the complex magnetic 
structure of actinide compound. For the purpose to improve the present theory from 
the quantitative viewpoint, it may be important to include also the deformation of local 
/-electron wavefunction, but we postpone such a study in the future. In this paper, 
we consider the splitting energy between r 8 orbitals, in order to take into account the 
tetragonality in the 115 compounds. 

2.2. r 8 model 

After tedious calculations, we obtain the Hamiltonian including T 8 orbitals in the form 
of HUHE?] 

H = H kin + H C ef + He, (1) 

where H kin denotes the kinetic term of 5/ electrons, Hqef is the CEF potential term, 
and He indicates the Coulomb interaction term among T 8 electrons. 

Concerning the kinetic term, one traditional way is to consider the hybridization 
process between / and conduction electrons. For 115 systems, from the band-structure 
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calculations (HEj , P electrons of Ga ion play an important role to form conduction band. 
However, based upon a picture that 5/ electrons are itinerant, we can consider another 
way to construct the kinetic term for / electrons within a tight-binding approximation. 
In actual compounds, there should occur several processes through other ligand ions, in 
addition to the direct process of / electrons. However, when we include correctly the 
local symmetry of relevant / orbitals, it is believed that we can grasp qualitative point 
of kinetic term of / electrons, even within the simple tight-binding approximation. This 
point will be discussed elsewhere in future. Here we do not consider the hybridization 
process. Then, the kinetic term is given in the tight-binding approximation for / 
electrons as 

#kin = 2j Kr'firafi+ar'a, (2) 
i,a,<7,T,T / 

where f\ Ta is the annihilation operator for an / electron with pseudospin a in the 
r-orbital at site i and t* T , is the /-electron hopping matrix element between r- and t'- 
orbitals along the a direction. Indices a and b denote the Tg and Tg orbitals, respectively. 
In the xy plane and along the z-axis, t* r , is given by 

e -t( 3/4 (3) 

-a/3/4 1/4 J' {6) 



for the x-direction, 



/ 3/4 VS/4\ 



^ v/3/4 1/4 J ' 



for the y direction, and 

for the z direction. Note that t—(3/7)(f fa), where (ffcr) is a Slater-Koster overlap 
integral between / orbitals through a bond. In the following, t is taken as an energy 
unity. Remark that we introduce another hopping amplitude along the z-axis. In 
AnTGas (An=U and Np), AnGa2 layer is sandwiched by two AnGa3 sheets, indicating 
that the hopping of / electron along the z-axis should be reduced from that in AnGa3. 
However, it is difficult to estimate the reduction quantitatively, since it is necessary 
to include correctly the hybridization with d electrons in transition metal ions and p 
electrons in Ga ions. Thus, in this paper, we consider the effective reduction by simply 
treating tl b as a parameter. For a practical purpose, we introduce a non-dimensional 
parameter t z =tl b /t. 

We point out that the hopping amplitudes among T 8 orbitals are just the same as 
those for the e g orbitals of 3d electrons (TJ |U] . Readers can intuitively understand this 
point from the shapes of Tg orbitals shown in Fig. 3(a), which are similar to e g orbitals. 
Mathematically, this is quite natural, if we recall the fact that T 8 is isomorphic to T 3 x T 6 , 
where T 3 indicates E representation for the orbital part and T 6 denotes the spin part. 
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Note, however, that the agreement between /- and (^-electron hopping amplitudes is due 
to the choice of special axis directions. If we take the (1,1,0) direction in an fee lattice, 
/-electron hopping amplitude becomes complex, leading to the stabilization of octupole 
ordering in Np0 2 [H2|- 

The CEF term is given by 

H CEF = -A^2( Pia - Pib )/2, (6) 

i 

where p\ TU = f\ Ta f\ T(J i Pir=J2a P\tm an d A is the level splitting, which expresses the effect 
of tetragonal CEF, as discussed above. We note that in actuality, A should be related 
to the value of t z , since both quantities depend on the lattice constant along the z axis. 
However, the relation between t z and A is out of the scope of this paper and thus, here 
we simply treat them as independent parameters. 
The Coulomb interaction term is expressed by 

i,r i 

+ J ^2 fiaafiba'fiaa'fiba + J' fljir'lfir^ , (7) 

i,o-,<r' i,r^T' 

where the Coulomb interactions U, U', J, and J' denote intra-orbital, inter-orbital, 
exchange, and pair-hopping interactions among T 8 electrons, respectively, expressed by 
using the Racah parameters [H31- Note that the relation U=U'+J+J' holds, ensuring 
rotational invariance in pseudo-orbital space for the interaction part. For rf-electron 
systems, one also has the relation J— J'. When the electronic wavefunction is real, this 
relation is easily demonstrated from the definition of the Coulomb integral. However, 
in the j-j coupling scheme the wavefunction is complex, and J is not equal to J' in 
general. For simplicity, we shall assume here that J=J', noting that essential results 
are not affected. Since double occupancy of the same orbital is suppressed owing to the 
large value of U, pair-hopping processes are irrelevant in the present case. 

3. Results 

Now let us show our calculated results based on the two-orbital model. First we discuss 
the electronic properties of if^m i n the non-interacting limit by focusing on the Fermi- 
surface structure, in order to gain an insight into the interpretation of dHvA experiments 
on Np-115 compounds. Next the magnetic properties are discussed. We attempt to 
understand the appearance of three kinds of magnetic phases based on the orbital-based 
scenario similar to that of manganites. 

3.1. Fermi- Surf ace Structure 

Let us consider the structure of Fermi-surface sheets of the r 8 tight-binding model -f/kin- 
In the following, we define n as the number of / electrons included in the T 8 state. 
Namely, the local /-electron number per actinide ion is n+2, by adding electrons. 
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(c) n=1 .5, tz=1 , A=0 (d) n=1 .5, tz=0.1 , A=-0.4 




Figure 4. Fermi-surface sheets of the Ts tight-binding model for (a) n=2, t z = l, A=0, 
(b) n=2, t z =0.1, A=-0.4, (c) n=1.5, t z =l, A=0, and (d) n=1.5, t z =0.1, A=-0.4. 
The bounding box indicates the first Brillouin zone for a simple cubic lattice. The T 
point is located at the center of the box, while the apices denote R points. 



As shown in Fig. 3(d), in our picture, there are two T 8 electrons in trivalent 
neptunium ion. Then, it is quite natural to consider first the case of n—2. The results 
for (t z , A) = (1,0) and (0.1,-0.4) are shown in Figs. 4(a) and (b), respectively. For the 
cubic case with t z —l and A=0, there exist two kinds of cube-like Fermi surfaces. One 
is centered at F point, while the center of another cube is R point. When we change 
t z and A to consider the tetragonality, cubes gradually change to cylinders. As shown 
in Fig. 4(b), we can clearly observe two kinds of cylindrical Fermi surfaces. Note that 
the Fermi-surface structure does not sensitively depend on the value of A, as long as 
the absolute value is not so large as |A|<0.5. It is important to have quasi orbital 
degeneracy, to reproduce a couple of cylindrical Fermi surfaces. 

It is interesting to recall the dHvA experiment on NpCoGas, which has revealed 
two kinds of cylindrical Fermi surface [52] • In the relativistic band-structure calculations 
for the paramagnetic phase, it has been difficult to understand the appearance of a 
couple of cylindrical Fermi surfaces jHH] ■ Note that the folding of the magnetic Brillouin 
zone cannot explain the discrepancy between the experimental results and the band- 
calculation ones for NpCoGas. A direct way to understand the dHvA experimental 
results is to improve the band-structure calculation method, but as shown in Fig. 4(b), 
even in a simple tight-binding model for itinerant 5/ electrons, we obtain two kinds of 
cylindrical Fermi surface due to the effect of r 8 orbital degree of freedom. We do not 
claim that the dHvA results can be understood only by the present simple calculations, 
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but a hint to understand the electronic structure of Np-115 is believed to be hidden in 
the construction of our model with r 8 orbital degree of freedom. 

We point out that in actual Np-115 compounds, the number of n is not definitely 
determined. Especially, due to the change of transition metal ion, the number of n seems 
to be changed effectively. In fact, in the relation between superconducting transition 
temperature T c and the ratio of lattice constants c/a, the curve of T c vs. c/a for solid 
solution (U,Np,Pu)CoGa 5 is similar to that for another solid solution Pu(Fe,Co,Ni)Ga 5 
[OH IBS]- Namely, in actinide 115 materials, the effect of the change of /-electron number 
due to the substitution of actinide ion seems to be similar to that of the substitution 
of transition metal ion. Thus, when we consider the change of the magnetic structure 
among NpTGas (T=Fe, Co, and Ni) based on the effective two-orbital model, there is 
no strong reason to fix firmly n—2. 

In Figs. 4(c) and (d), we show the Fermi-surface sheets for n=1.5, in which the 
number of / electrons is slightly decreased. We can observe that the Fermi-surface 
structure is not changed so much, although the cylinder centered at the F point becomes 
slender. Thus, as long as we are based on the simple tight-binding model, it is concluded 
that the structure of the Fermi-surface sheets is not changed sensitively by the number 
of n around at n=2. This result suggests that it is possible to change the value of n to 
consider the magnteic structure of Np-115 compounds. 

3.2. Magnetic structure 

Now we consider the magnetic properties of the r 8 model. Among several methods 
to analyze the microscopic model, in this paper we resort to an exact diagonalization 
technique on a 2x2x2 lattice. Although there is a demerit that it is difficult to enlarge 
the system size, we take a clear advantage that it is possible to deduce the magnetic 
structure by including the effect of electron correlation. In order to discuss the ground- 
state properties, it is useful to measure the spin and orbital correlations, which are, 
respectively, defined by 



with T*=Y^ a (n\MT— 7"H&(t)/2. Here N is the number of sites. 

First let us consider the case of n=2, which is corresponding to the trivalent Np 
ion. In Fig. 5(a), we show the spin correlation as a function of A for U'—5, J—0, 
and t z =l. At t z =l and A=0 (cubic system), local triplet composed of a couple of / 
electrons is formed at each site and the G-type AF structure is stabilized due to the 
so-called superexchange interaction. As easily observed, even when A is introduced as 
the tetragonal CEF effect, the G-AF structure remains robust for |A| < 1. When |A| is 
larger than unity, two electrons simultaneously occupy the lower orbital, leading to the 




(8) 



with crf=5^ r (n iTl — n iTi )/2, and 




(9) 




Figure 5. Spin correlation for n=2, {7'=5, and J=0 as a function of (a) A for t z =l 
and of (b) t z for A=0. 



non- magnetic state composed of local r 1; irrelevant to the present study to consider the 
magnetic phase. In Fig. 5(b), the spin correlation is shown as a function of t z for A=0. 
Again, the G-type AF structure is stabilized, but we find that the spin correlation of 
qr=(7r, 7T, 0) comes to be equivalent to that of q=(ir, n, ir) with the decrease of t z , since 
the AF structure is stabilized in each xy plane due to superexchange interaction and 
the planes are decoupled for small t z . 

At the first glance, it seems difficult to understand the variety of magnetic phases 
observed in NpTGas, when we consider only the trivalent Np ion. However, there is no 
a priori reason to fix the valence as Np 3+ , as mentioned in the previous subsection. In 
NpTGas, (^-electron band originating from transition metal ions may significantly affect 
the valence of Np ion. In addition, we also stress that the actual compounds exhibit AF 
metallic behavior. In the band-structure calculation, the average number of / electrons 
at Np ion is easily decreased from four. Thus, we treat the local /-electron number as 
a parameter. 

We may consider another reason to decrease effectively the number of / electron 
from n=2 in NpGa3. In the present two-orbital model, the G-AF structure is robust, 
which is natural from the theoretical viewpoint within the model. However, in the 
experimental result on NpGa3, the low-temperature ground state is ferromagnetic, 
although the AF phase has been observed around at T~60K [66J. In order to understand 
the occurrence of the FM phase in the two-orbital model, it is necessary to inject some 
amount of "hole" in the AF phase, since the double-exchange mechanism works to 
maximize the kinetic motion of electrons, as will be discussed later in the main text. It 
is difficult to determine the amount of doped holes to obtain the FM phase, but at least 
qualitatively, the effective decrease of n seems to be physically meaningful in NpGa3 as 
well as in NpTGas. 

Now we consider the case of n=1.5. As discussed in the previous subsection, the 
Fermi-surface structure does not change so much from that of the case of n—2, when 
| A| is not so large. In Fig. 6(a), we show the ground-state phase diagram in the (U', J) 
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Figure 6. (a) Ground-state phase diagram in the (£/', J) plane for n=1.5, t z =l, and 
A=0. (b) Spin structures in the G-type AF, C-type AF, A-type AF, and FM phases. 

plane at n—1.5 for the cubic case with t z =l and A=0. At J=0, a G-type AF structure 
is stabilized due to superexchange interaction in the same way as the case of n=2. 
However, the G-AF structure is immediately changed to a C-AF structure only by a 
small value of the Hund's rule interaction. With increasing J, the magnetic phase 
changes in the order of G-AF, C-AF, A-AF, and FM phases, except for the C-AF phase 
in the large J region. The spin structure of each phase is shown in Fig. 6(b). This 
result is quite natural, since we are now considering the magnetic structure based on 
the two-orbital model, in which the FM tendency is due to the optimization of kinetic 
motion of electrons. 

In order to understand the underlying mechanism of the appearance of various 
magnetic structures in addition to the G-AF structure for n—1.5, we consider the case 
of J—0. In Fig. 7(a), we show typical results of the spin correlation as a function of A for 
U'—5, J—0, and t z —l. At A=0, the dominant component of S(q) appears at q={ji, ir, n), 
indicating the G-type AF structure. On the other hand, in the region of positive small 
A, the spin correlation of q=(ft, vr,0) turns to be dominant, indicating the C-type AF 
structure. This phase is defined as C-AF(I). Concerning the orbital correlation, we find 
no remarkable structure in the C-AF(I) phase, as shown in Fig. 7(b). Note that in 
the G-AF phase with positive A, the orbital correlation of qr=(0,0, 0) has a significant 
value in comparison with other regions, since electrons tend to doubly occupy lower Tg 
orbitals. This G-AF phase is called G-AF(I). On the other hand, in the G-AF phase 
with negative A, there is no indication of the development of the orbital correlation of 
q—(0, 0,0), since the lower Tg orbitals are singly occupied. We label this G-AF phase 
as G-AF (II) to distinguish it from G-AF (I). It is interesting to see different types of the 
G-AF phases due to the change of the orbital state. 
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Figure 7. (a) Spin and (b) orbital correlations as a function of A for n=1.5, U'=5 
J=0, and t z — l. (c) Schematic view of the electron configuration in the C-AF phase 
in the region near t z =l. (d) Ground-state phase diagram of the magnetic structure in 
the (A,t z ) plane. 

Let us now discuss the reason of the appearance of the C-AF (I) phase in the region 
of positive small A. For A>0, the energy level of Tg is lower than that of Tg by definition. 
Then, the Tg orbital is first occupied by one electron at each site. The rest of electrons 
are accommodated in Tg orbitals to avoid the effect of intra-orbital Coulomb interaction 
when A is not large so much. The electrons in Tg orbitals hop to the Tg orbital at the 
nearest-neighbor site, since the T b 8 has higher energy. Since the electrons in Tg orbitals 
move only in the xy plane, the AF structure is stabilized in each xy plane independently 
due to superexchange interaction, as schematically shown in Fig. 7(c). On the other 
hand, the electrons in Tg orbitals can move itinerantly within the network composed of 
Tg orbitals due to the existence of holes in Tg orbitals. Since the hopping amplitude 
along the z axis is dominant in the region near t z =l, the motion of Tg electrons causes 
an FM spin arrangement of Tg electrons along the z axis to gain kinetic energy. Namely, 
even though localized t2 g spins do not exist in the present case, the so-called double- 
exchange mechanism is always active in the two-orbital model. Note that the C-AF 
structure disappears for large A, since all the electrons tend to occupy orbitals and 
the double-exchange mechanism is no longer active in that situation. 
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Figure 8. Ground-state phase diagram of the magnetic structure in the (A, t z ) plane 
for n=1.5, U'=5, (a) J=0.5 and (b) J=4. 

After calculations of the spin and orbital correlations for several parameter sets, 
we obtain the ground-state phase diagram in the (A,t z ) plane, as shown in Fig. 7(d). 
It is found that the C-AF(I) phase discussed above extends to a wide region of t z <\ in 
the phase diagram. Note that for small values of t z , another C-AF phase, which we call 
C-AF(II), appears in the region of negative small A, in which the spin correlations of 
(ir, 0, 7r) and (0, ir, ir) are dominant. For small t z and large negative A, there appears yet 
another G-AF phase, called G-AF(III), due to the same reason as that for the occurrence 
of G-AF(I) phase for positive A, but the occupied orbital is changed from Tg to r|. In 
any case, for J=0, we can observe several kinds of C- and G-AF phases, but A-AF phase 
does not occur. 

Although we increase the value of J as J— 0.5, no new phases appear in the phase 
diagram for n—1.5, as shown in Fig. 8(a). There are three phases, but they are two 
C-AF and one G-AF states. As labelled explicitly in the phase diagrams, C-AF (I), 
C-AF(II), and G-AF(I) are the same as those in the phase diagram of Fig. 7(d). Due 
to the effect of J, G-AF (II) and G-AF (III) disappear, since the number of FM bonds 
should be increased to gain the kinetic energy. As shown in Fig. 8(b), when we further 
increase the value of J as J=4, the G-AF phase completely disappears and instead, we 
observe the A-AF phase sandwiched by two C-AF phases. As described above, due to 
the double-exchange mechanism in the two-orbital model, the A-AF phase is considered 
to appear, when J is increased. 

4. Discussion and Summary 

In this paper, we have proposed an effective model with active orbital degree of freedom 
to understand the magnetic structure of neptunium 115 compounds from the microscopic 
viewpoint. By analyzing the model with the use of the exact diagonalization technique, 
we have obtained the ground-state phase diagrams including three kinds of AF phases 
corresponding to NpTGa 5 . 
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In the experiments for NpTGas, C-, A-, and G-AF magnetic phases have been 
found in NpFeGas, NpCoGas, and NpNiGas. Here we have a question: What is a key 
parameter to understand the change of the magnetic structure? In the case of UTGa5, 
it has been claimed that the level splitting A is important to explain the difference in 
magnetic structure as well as the magnetic anisotropy for a fixed value of n=l [57J. 
Roughly speaking, A is positive for T=Fe, small positive for T=Co, and negative for 
T=Ni. Among UTGa 5 with T=Ni, Pd, and Pt, when we assume that the absolute 
value of A is increased in the order of Ni, Pd, and Pt, it is possible to understand 
qualitatively the change in the magnetic anisotropy, in addition to the change in the 
magnetic structure of G-AF for T=Ni and A-AF for T=Pd and Pt. It has been found 
that the value of t z is not so crucial to explain qualitatively the magnetic properties of 
U-115 based on the two-orbital model for n—1. 

For n=2, as emphasized in the previous section, we always obtain the G-AF phase. 
However, for n=1.5, we have observed three kinds of AF magnetic structure in the 
phase diagrams. Let us summarize the change in the magnetic structure for a fixed 
value of t z =0.8. Note that this value is larger than t z —0.1, which we have considered 
to reproduce two kinds of cylindrical Fermi-surface sheets of Np-115. However, in the 
small-sized cluster calculations, it is difficult to compare directly with the values in the 
thermodynamic limit. Thus, we do not discuss further the quantitative point on the 
values of t z here. As shown in Fig. 7(d), for J=0 and ^=0.8, we see the change in the 
magnetic structure as G-AF (A<0), C-AF(0<A<0.4), and G-AF (A>0.4). For J=0.5 
and t z =0.S, as shown in Fig. 8(a), the C-AF phases are always observed, but they have 
different orbital structures. Finally, for J=4 and t z =0.8, we observe C-AF (A<-0.15), 
A-AF(-0.15<A<0.3), and C-AF (A>0.3), as shown in Fig. 8(b). 

In order to understand the appearance of three types of the AF phases, we may 
consider an explanation due to the combination of the changes in A and n. For instance, 
by assuming that J=4 for NpTGas and the change in A for NpTGas is just the same 
as that for UTGas, we consider that n^2 with A<0 for T=Ni, n~1.5 with A^O for 
T=Co, and n~1.5 with A> for T=Fe. Then, it seems to be possible to relate our 
theoretical AF phases with the experimental observations in NpTGa 5 . However, it is 
difficult to claim that the above parameter assignment for three Np-115 materials is the 
best explanation for the magnetic structure of Np-115, since in actual compounds, there 
are other important ingredients which have not been included in the present model. 
For instance, we have never discussed the direction of the magnetic moment of Np ion. 
Especially, the canted AF structure cannot be considered at all for the G-AF phase of 
NpNiGas. Thus, we need to recognize some distance between the actual magnetic states 
and the theoretically obtained phases. Our theory should be improved by taking into 
account other ingredients of 115 structure. 

In summary, we have analyzed the orbital degenerate model appropriate for 
NpTGa 5 by using the exact diagonalization technique. Our phase diagram includes C-, 
A-, and G-AF phases. We have proposed the manganite-like scenario to understand the 
appearance of three kinds of AF spin structure in Np-115. Namely, the double-exchange 
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mechanism works also in some actinide compounds based on the model with active 
orbital degree of freedom. We believe that the present model can grasp some important 
points of the actinide compound by regarding it as charge-spin-orbital complex. 
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